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Abstract 

We show that the manifestation of quantum interférence in graphene is very différent from that in 
conventional two-dimensional Systems. Due to the chiral nature of charge carriers, it is sensitive not 
only to inelastic, phase-breaking scattering, but also to a number of elastic scattering processes. We 
study weak localization in différent samples and at différent carrier densities, including the Dirac 
région, and find the characteristic rates that détermine it. We show how the shape and quality of 
graphene flakes affect the values of the elastic and inelastic rates and discuss their physical origin. 
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The quantum correction to the conductivity of two-dimensional Systems due to électron 
interférence has been studied for more than twenty years [U This phenomenon of weak 
localization (WL) has become a tool to détermine the processes responsible for électron 
dephasing due to inelastic électron scattering or scattering by magnetic impurities [21 E]. 
In this well-established field of research it comes as a surprise to discover that in a new 
two-dimensional System, graphene WL does not follow the standard convention that it 
is only controlled by inelastic and spin-flip processes. First attempts to measure WL in 
graphene have produced contradictory results that tentatively point towards this unusual 
behavior El [7j. Measurements on graphene flakes fabricated by mechanical exfoliation 
[5] have shown that in the majority of samples WL is totally suppressed. In contrast, in a 
sample fabricated by an alternative, epitaxial method, WL has been distinctly seen, albeit 
at a single (high) carrier density f7\. 

The theory of WL in graphene [8] predicts a remarkable feature: it should be sensitive not 
only to inelastic, phase breaking processes, but also to différent elastic scattering mechanisms 
[SI m [TU]. The reason for this is that charge carriers in graphene are chiral, that is, they 
have an additional quantum number (pseudospin) [TIJ. Elastic scattering that breaks the 
chirality will destroy the interférence within each of the two graphene valleys in /c-space. 
Such defects, characterised by the scattering rate r~^, include surface ripples, dislocations 
and atomically sharp defects [5l jTU] . Intra-valley WL can also be destroyed by anisotropy 
of the Fermi surface in fc-space, so called 'trigonal warping' [S], characterised by the rate 
T^^. There is one elastic process, however, which acts to restore the suppressed interférence. 
This is mier- Valley scattering, which occurs at a rate Tj"^ on defects with size of the order 
of the lattice spacing a. As the two valleys have opposite chirality and warping, inter- 
valley scattering is expected to negate the chirality breaking and warping effects by allowing 
interférence of carriers from différent valleys. 

In this work we aim to examine what factors are responsible for the manifestation of WL in 
graphene fabricated by mechanical exfoliation [1] . We study the magnetoconductivity (MC) 
in perpendicular field of several samples with différent quality and dimensions, with the aim 
to control the relation between the scattering rates of carriers. Thèse studies are performed 
at différent carrier densities controlled by a gâte voltage Vg, which include densities around 
the Dirac point at = 0. This point is spécial as about it there is a change of the type of 
carrier from électrons to holes and therefore the net carrier density is zéro. The conductivity. 
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however, is seen to remain at a finite value (Tmin ~ /h and not drop to zéro [12 
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FIG. 1: Resistivity of graphene flakes as a function oi Vg ai T =0.25 K. The mobilities (in 
(zm^\-^s-^) of the samples outside the Dirac région: 5100 (D), 7500 (Fl), 10000 (F2) and 8000 
(B). The insets show the first quantum Hall plateau, where filling factor v = nh/4eB. The right 
panel shows SEM images of the samples, where the positions of the contacts are shown as outlines. 
The diagram at the foot of this panel shows the graphene sample on n+Si substrate (purple), cov- 
ered by 300 nm Si02 (blue) and contacted by Au/Cr (yellow). Control of the carrier density n is 
achieved by Vg. 



Figure 1 shows the resistivity as a function of Vg of four samples with différent shapes 
and mobilities: D, Fl, F2 and B, each with a typical peak around the Dirac point. Sample 
D is a square flake; Fl and F2 are rectangular with similar width to, but length larger than 
D; B is a narrow strip of similar length to Fl and F2 but with much smaller width, Fig. 1 
(right). Insets in Fig. 1 demonstrate measurements of the first quantum Hall plateau, which 
shows a half-integer step (0.5 x 4e^//i) - a clear indication that the samples are single-layer 
graphene [12]. We wanted to see what différence the shape of the samples will make to the 
WL - e.g., the narrowest sample B is expected to have the largest scattering rate r~^, as the 
edges could produce strong inter-valley scattering. To understand the relation between the 
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scattering and the détails of the graphene surface, the electrical measurements have been 
complemented by atomic force microscope (AFM) imaging of the sample topography, which 
have shown the présence of ripples, and additional rapid folds (ridges) across sample B (see 
Supplementary Material below). 
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FIG. 2: (color online) Magnetoconductivity observed in graphene flakes. (a) Dirac région of sample 
1^1 ^ 1 V, n < 7 X 10^°cm~^; (b) sample D, ~ 14 V, n ~ lO^^cm"^ (the legends of (a) and 
(b) are the same); (c) samples Fl and F2 at T = 1 K, ~ 10 V, n ~ 7 x 10^^ cm~^; (d) sample B, 
Vg ~ 11 V, n ~ 8 X 10^^ cm^^. Lines are best fits to Eq. [ï| 

In order to study the conductivity correction caused by WL we must first account for 
the reproducible conductance fluctuations clearly seen in Fig. 1. They are also présent as a 
function of magnetic field B and are caused by the fact that the graphene flakes are small 
(comparable to the dephasing length L^), so that the average effect of WL is of the same 
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order as the fluctuations which have the same, quantum interférence origin [2]. We use tlie 
procédure developed in [13] in tlie study of WL in bilayer grapfiene: RiVg) is measured at 
différent B and the results are averaged over a range AVg = 2 V shown in Fig. 1 by boxes 
which contain a large number of fluctuations (see Supplementary Material below). Examples 
of the obtained average MC, Aa{B) = {a{Vg,B) — cr(yg,0)) AVg, are shown in Fig. 2 for 
différent samples. For the analysis of the results we use the theory [8.j which predicts that 
the MC is controlled by several scattering rates, both inelastic (t^^) and elastic (t-~^, r~^, 

Here F(z) = Inz + tjj(0.5 + z^^), tp^x) is the digamma function, r^^ = AeDB/h, r^^ is 
the phase-breaking rate and t^^ = t^^ + t^^. (The theory assumes that the momentum 
relaxation rate r^^ is the highest in the System and comes from the Coulomb charges in 
the Si02 substrate and not atomically sharp defects, so that it does not affect the carrier 
chirality.) The first term in Eq. (1) is responsible for weak localization, while the terms 
with négative sign can resuit in 'anti-localization'. 

We have found that among the différent possible relations between the scattering rates, 
the foUowing holds in ail studied samples: the intra-valley WL is strongly suppressed due 
to a large rate t^^, which approaches Tp^; however, WL is clearly seen in ail régions of the 
carrier density, due to significant intra-valley scattering, r^~^ ~ t^^. At the same time, the 
shape of the MC curves can be very différent as it is controlled by the interplay between 
ail scattering rates involved, Fig. 2. Comparing two régions of carrier densities for square 
sample D (Dirac région (a) and électron région (b)) one can see that in (a) the curves have a 
much stronger downturn, indicating greater importance of the third ('anti-localising') term 
in Eq. [î] due to smaller rate r~^. For two geometrically similar samples Fl and F2 in 
Fig. 2(c), it is seen that sample F2 (with largest mobility) has a more rapid increase of 
the conductivity in smaller field (due to smaller r^^) and more rapid downward turn of the 
curves at larger fields (due to smaller r^"^). For the narrow sample B, Fig. 2(d), the MC 
curves do not turn down at ail, indicating a very fast inter-valley rate r^~^ and therefore 
unimportance of ail terms in Eq. 1 apart from the first. 
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FIG. 3: (color online) Characteristic lengths responsible for weak localization; dotted lines are 
guides to the eye. Sample D: (a) the Dirac région {n < 7 x W^^ cm^^) and (b) électron région 
(n ~ 10^^ cm~^); (c) Phase-breaking rate t^^ = D/L'^ as a function of T for différent n. Inset to 
(a) illustrâtes the saturation of ai low T due to the sample size. Inset to (b) shows the scattering 
process behind the length Lj. Sample F2: (d) Température dependence of the characteristic lengths 
in the électron région (n ~ 10^^ cm^^). 



Figure 3 shows the température dependence of the characteristic lengths found from the 
analysis of the MC by the best fit with Eq. [î] Figure 3(a,b) compares the results for the 
Dirac and électron régions for sample D, where it is seen that is température dépendent 
at high T (> 3 K) but saturâtes at a value L^^* at low températures. Figure 3(b,d) compares 
the results for samples D and F2, at close values of carrier densities outside the Dirac région. 
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Sample F2 is about three times longer, and one can see that L^^^ is larger in the longer sample 
F2. (In sample B, L^^^ has also been found to be larger than in sample D.) This clearly 
implies that the reason for the saturation is a limitation imposed by the sample size, and 
not by scattering from a small, uncontrolled number of magnetic impurities |3]. 

In the Dirac région L^^* has been found to have a smaller value than at larger carrier 
densities. (It is interesting to note that the narrow sample, B, shows the biggest decrease 
of L^J^^ in the Dirac région, while the square sample, D, the smallest.) This decrease can be 
related to the inhomogeneity of the sample at low carrier densities. It can resuit in formation 
of electron-hole puddles, so that at = there are equal densities of électrons and holes 
and not zéro density of each type of carrier. Inhomogeneity can modify the geometry of 
conducting paths and decrease the effective dimensions of the sample, resulting in a smaller 
value of L^^^. 

The température dependence of contains information about the inelastic mechanism 
responsible for the dephasing of charge carriers. There are suggestions that electron-electron 
(e-e) interaction, the main mechanism of dephasing at low T, is différent in graphene com- 
pared with other Systems [II]. To examine this we have analyzed the T— dependence of the 
dephasing rates found from analysis of the WL. Figure 3(c) shows the phase-breaking rate in 
différent density régions of sample D. (To find t^^ we use the relation = (Dr^Y^"^, where 
the diffusion coefficient D = vpl/'i is determined from the mean free path / = h/2e^kFP- 
For the Dirac région, where the puddles can be formed, the value of the Fermi wavenumber 
kp (inside the puddle) is simply estimated at the boundary of the région where the MC is 
studied, = 1 V, Fig. 1.) Our results show that électron dephasing rate obeys the usual, 
linear T— dependence for e-e scattering in the 'dirty limit', Tr^ < 1 [Ij: r^^ = pkBTlïig/hg, 
where g = ah/e^. (In our samples the parameter Txp varies from 0.002 to 0.4 in the studied 
température range 0.25-25 K.) The empirical coefficient (3 is found to be between 1 and 
2 in ail studied régions, Fig. 3(c), and ail samples. Therefore, we can conclude that while 
électron interférence in graphene is significantly différent from other Systems, e-e interaction 
does not show unconventional behavior. 

In addition to the analysis of the WL, we have also analyzed conductance fluctuations as 
a function of B and Vg using standard relations in terms of L,^ [2]. This analysis has given 
values of close to those found from the analysis of WL (see Supplementary Material 
below) . 
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Now let us discuss the behavior of the elastic, inter-valley length Li which we have verified 
to be essentially T— independent in ail samples. In samples D, Fl and F2 the found Lj is 
comparable to the width of the samples (approximately half the width). This means that, 
indeed, the sample edges make significant contributions to inter-valley scattering. This is 
consistent with the fact that the narrowest sample B has shown the smallest value of Lj. 
However, the value of Li for sample B is about three times smaller than the sample width. 
This can be due to the présence of rapid ridges of height ~ 1.5 nm observed in this sample 
by AFM, Fig. 4(a). They are separated by a distance smaller than the sample width and can 
be another source of inter-valley scattering. This suggests that the inter-valley scattering is 
controUed not only by the edges but also by the defects in the inner part of the sample. 

We have found that the intra-valley scattering length L^, is much smaller than the inter- 
valley length Li, Fig. 3, and approaches the mean free path. There are several possible 
mechanisms that can be responsible for the observed large intra-valley scattering rate and 
resulting strong suppression of WL in one valley (see Supplementary Material below). Scat- 
tering by atomically-sharp defects cannot explain this: such scattering is also a source of 
strong inter-valley scattering, so that Li and would be comparable if this mechanism was 
dominant. The smaller value of in experiment must therefore be due to an additional 
scattering rate which affects L^, but not Lj: from warping |8|, or from the defects of the 
crystal structure that are large on the atomic scale O HU]. Estimation of the expected 
using theory [S] gives a value of < 0.3 ps^^ which is much smaller than the expérimental 
10 ps ^. Therefore, the reason for small L^, could lie in the defects of the crystal 
structure of graphene flakes: ripples and dislocations. (The strain in the lattice induced by 
such defects acts as a source of effective magnetic field that can destroy WL.) As the dislo- 
cation core is also a source of inter-valley scattering, their séparation can be estimated from 
the known value of Lj. This value is much larger than the dislocation séparation ^ ~50nm 
required to explain the small value of L^,. (This value of ^ is obtained using the relation 
r^^ = vp/kp^'^ [lÏÏ]-) The effect of ripples on the graphene surface is also negligible in 
our samples, if we use the estimation of this effect from [5]. The roughness of our samples 
found from AFM measurements, Fig. 4(b,c), is /i ^ 0.3 nm and d ^ 10 nm (in agreement 
with [15j), which gives for a typical L^p ^ 1 /im an effective magnetic field B^s ~ 1 niT. 
This is a small correction to the real fields used in experiment, Fig. 2. There is one more 
mechanism that can introduce an asymmetry in the crystal structure and hence break the 
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chirality of carriers: a potential gradient coming from charged impurities in the substrate. 
Our estimation of its effect using the approach of [TU] has also given a negligible resuit (see 
Supplementary Material below). Therefore, the obtained values of the scattering rate 
are much higher than those predicted by existing models, and more detailed théories of the 
mechanisms of intra-valley suppression of WL in graphene are required. 
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FIG. 4: Atomic force microscope images of graphene surfaces. Below the image of sample B (a) is 
the surface profile averaged over the width of the sample. Below the images of the topography of 
samples D (b) and F2 (c) are the corresponding autocorrélation functions of the surface roughness. 

In summary, we have shown that the weak localization correction in graphene exists at 
ail studied carrier densities, including the Dirac région. Its manifestation is determined by 
the interplay of inelastic and elastic scattering mechanisms, which makes WL a sensitive 
tool to detect the defects responsible for inter-valley scattering and chirality breaking. We 
show that, in spite of a strong intra-valley suppression of WL, the quantum interférence 
correction to the conductivity is clearly seen due to significant inter-valley scattering. Total 
suppression of WL is only possible in experiments where inter-valley scattering is negligible, 



i.e. in very large samples without sharp defects in the bulk. 
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Samples 

Samples were manufacturée! using the method of mechanical exfoliation of highly-oriented 
pyrolytic graphite devised in [3], on a n"'"Si/Si02 substrate with oxide layer of thickness 
t = 300 nm. Lithographically defined Au/Cr contacts were subsequently made to each flake. 
Résistance measurements were carried out in the température range from 0.25 to 25 K using a 
standard lock-in technique with 1 nA driving current. Samples B, D, Fl are two-terminal and 
F2 is four-terminal (the additional contacts were used to account for the contact résistance). 
The concentration of carriers (électrons n and holes p) in graphene is determined by the 
capacitance between the graphene and n+Si substrate: e{p — n) = (eeo/t)V^. There was a 
small unintentional doping of the samples leading to a shift in gâte voltage (~ 5V) of the 
position of the résistance peak with respect to = 0, which has been accounted for in the 
main text. The graphene-Au/Cr contact résistance has been found from the déviation of 
the height of the quantum Hall plateau from the expected value of 2e^ /h (see insets to Fig. 
1 of main text). The values of the contact résistance for samples Fl and D are about ~100 Vl 
and ~600r2 for sample B. 

Averaging procédure and analysis of magnetoconductance 

A method of effective averaging is important in small-sized samples to remove the in- 
fluence of mesoscopic fluctuations, as without it one can get contradictory results for the 
magnetoconductance (MC). (If we attempt to measure Aij{B) at différent V^, the character 
of the MC dépends on the spécifie point in Vg at which it is measured). Figure [î] shows 
how the averaging is performed. For each température the conductivity of the sample as a 
function of the gâte voltage is first measured across a 2 V range at incrémental values of the 
magnetic field. Then the curve at zéro magnetic field is subtracted from each curve and the 
resulting différence is averaged across the 2 V gâte voltage range. One can see from Fig. 1 
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FIG. 1: Illustration of the averaging procédure of the magnetoconductance of sample Fl in two 
density régions at two températures (only a fraction of AVg is shown here): (a) Dirac région, (b) 
électron région. Dotted lines show repeated sweeps at B = 0. 

the average increase of (Acr)AVg with magnetic field. Thèse averaged values of the MC are 
plotted as a function of B in Figure 2 of the main text. 

The perturbation theory of weak locahzation (WL) is apphcable at kpl ^ ^ (a diffusive 
métal). In our samples kpl, found from the conductivity a = 2e'^{kFl)/h, varies in the range 
3-30, with the smallest values in the Dirac région: 4, 3, 8 and 6 for samples D, Fl, F2 and B, 
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respectively. Another limitation for the application of the diffusive theory of WL is i? < Btr, 
where the 'transport' magnetic field is found from the condition Lb = {h/eBY^"^ ^ /. This 
limits the range of magnetic fields where we perform the analysis to -B ^ 100 mT. 

For the narrowest sample B, the dephasing length is larger than its width and therefore 
the ID theory of WL [8] should be used in the analysis of its MC in small fields. However, 
at fields where Lb <W {W is the width of the sample) i.e. at S > 7 mT, the 2D theory 
becomes applicable. As the bulk of the data is obtained in this range of the field, we have 
used 2D theory (Eq. 1 in main text) to analyse the MC. 

Comparison of characteristic lengths and times 

Figure [2] shows for samples Fl, F2 and D a comparison of the length L^, with length Lj, 
as well as the values of the corresponding times and (using = {Dt^Y^"^) for différent 
carrier densities. We emphasise that in the analysis of the MC the value of L^, is closely 
linked to that of Lj. In Eq. 1 the second and third terms have the same sign, therefore by a 
slight increase of one of them and a corresponding decrease of the other, one can get a similar 
agreement with experiment. Figure [2] shows not only the values found from the best fit (the 
higher 5-region being most sensitive to thèse two parameters) but also the synchronous 
variation allowed in thèse values while retaining a good fit, indicated by arrows. In spite of 
the variations, there are several trends seen in the figure. First, the value of Lj is always 
significantly larger than and somewhat larger in the better quality sample F2. Second, 
there is a decrease of Lj with increasing carrier density, although its value is smaller in the 
Dirac région. Finally, there is a decrease of when the carrier density is increased above 
the Dirac région. The dashed curves in Fig. 2 indicate the expected decrease of Lj and if 
the scattering rate is proportional to the density of states, which increases linearly with the 
Fermi energy e^? oc Vg^"^ . 

Estimations of the effects suppressing WL in a single valley 

Trigonal warping 

According to [8] the breaking of the time-reversal symmetry in one valley can occur due to 
the suppression of backscattering by the trigonal warping of the Fermi surface. The trigonal 
warping rate is 

= 2Tp{fiel/hvlf , 

where Tp is the momentum relaxation time, vp ~ 10^ ms~^ is the Fermi velocity and n is the 
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FIG. 2: Comparison of characteristic lengths and times for samples Fl, F2 and D at différent 
carrier densities. 

structural parameter equal to /i = ^^a^ j'èf?. Hère 70 ~ 3eV is the nearest-neighbour hop- 
ping energy and a ^ 0.26 nm is the lattice constant in graphene. For the typical parameters 
in our samples we obtain ^ 0.001 ps^^ for the Dirac région [ep ~ 30meV, Tp ^ 0.1 ps) 
and T^^ ^ 0.3 ps"^ for the highest measured concentration {ep ~ 130 meV, Tp ^ 0.05 ps). 
Trigonal warping of the Fermi surface is therefore a very weak effect compared to other intra- 
valley scattering mechanisms and cannot be the main reason of the strong chirality-breaking 
observed in our experiments (r^^ ~ Tp~^)- 



Dislocations 
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Another possible mechanism of chirality breaking in the graphene sheet is dislocations 
in the honeycomb lattice [TU]. If the trajectory of a quasiparticle goes near the core of 
a dislocation it leads to a change of the phase due to the induced strain. For randomly 
distributed dislocations the scattering rate related to this mechanism is 



Vf 



sauge ^^^2 ' 

where ^ is the average distance between dislocations [TU]. In order to obtain the experi- 
mentally found chirality-breaking rate t^^ ^ 10 - 20ps~^ the distance ^ should be about 
15 — 50 nm. However, the cores of the dislocations should also cause inter-valley scatter- 
ing, which is why this estimation is in contradiction with the relatively large value of the 
inter-valley scattering length (Lj ^ 1/im) observed experimentally. 

Ripples 

As proposed in [S], ripples in the graphene layer on a silica substrate can lead to sup- 
pression of weak localization because of the effective magnetic field generated by strain of 
the interatomic bonds. The vector potential corresponding to a single ripple with diameter 
d and height /i is [5J: 

_ 70 IV/^I' ^ 

evF ' 

where Vh ^ h/d. The flux through one ripple is $ = ^ A ■ dl ~ Ad and $ = J B ■ dS ~ Bd^, 
therefore the magnetic field associated with one ripple is 

d evF d^ 

Since the curvature vector of a ripple is random, the resulting magnetic field through the area 
limited by the dephasing length and containing ^ L'^/d'^ ripples should be averaged 
as follows: 

B lo fh^^ 

-Deff = 



/N cvfL^ \d^ 

The roughness of the graphene sheet found from AFM measurements is about 0.3 nm and 
the size of the features is about 10 nm. This gives a value for the magnetic field associated 
with one ripple B ~ 0.1 T. For our typical value of L^", ~ 1 fim the effective magnetic field is 
then Bcs ~ 1 niT. Since suppression of the quantum interférence requires a magnetic field 
Bcs > -Btr ~ 0.1 T, the estimated value is too small to destroy the localization effect. The 
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random magnetic field can only introduce an uncertainty in the value of B, Fig. 2 of the 
main text, comparable to the accuracy to which the field is set by the power supply. 

Potential gradients 

The last mechanism which can produce the breaking of time-reversal symmetry is a 
gradient of potential coming from the charged impurities in the substrate. A potential 
gradient leads to a distortion of the dispersion curve of a single valley and hence breaks the 
Valley symmetry. As shown in [10] the resulting scattering rate can be estimated as 

In order to get r^j-l^^ ~ one should have kpa ~ 1. This corresponds to the carrier density 
n = kp/n ^ 5 - lO^^cm"^, which is two orders of magnitude higher than the densities studied 
in the experiment. 

We conclude from thèse calculations that ail existing estimations for the chirality-breaking 
scattering rates are not sufficient to explain our expérimental results. 

Scanning probe microscopy studies 

The atomic force microscope used in this work was an Ntegra Aura from NT-MDT. 
We used non-contact tips NSGOl with résonance of 150 Hz at an amplitude of < 40 nm. 
To obtain high resolution in the xy-plane 'diamond-like carbon' coated tips with curvature 
radius 1 — 3 nm were used; tip convolution therefore limited feature resolution to this scale. 
To remove the infiuence of the water layer présent on the silica substrate ail measurements 
were performed in an atmosphère of dry nitrogen at 3 mbar, giving a tip résonance quality 
factor Q ^ 1000. The noise in the 2;-scale (height) is of the order 0.02 nm measured on pure 
graphite and silica with the AFM operating with acoustic and vibrational isolation. 

We found that the surfaces of the silica and the graphene after lithographie processing 
were covered in droplets of PMMA with height ~ 2 nm, similar to the findings of [15] . They 
reduced the image quality and also made détermination of the step edge between graphene 
and silica difficult. To obtain the scans of clean graphene shown in Fig. 4 of the main text 
we mechanically cleaned the surfaces. Figure [3]^a) shows a phase contrast image of sample 
F2 where both the PMMA droplets and a cleaned area are seen. 

To understand the extent to which the PMMA droplets exist under the fiake (due to the 
lithographie process of depositing location markers prior to the déposition of the graphene 
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FIG. 3: AFM measurements of sample F2. (a) Phase contrast image where the PMMA droplets 
on the left and cleaned area on the right can be seen, with two régions in the clean area highlighted 
by boxes, (b) Magnified topographie image of box '(b)' showing the torn edge of the graphene 
flake. (c) Magnified phase contrast image of box '(c)' with silica (top) and graphene (bottom). (d) 
Autocorrélation analysis of the roughness in the boxes highlighted in (c), with insets of silica (top) 
and graphene topography (bottom). Scan (a) is 3 fim size and the phase change at the graphene- 
silica boundary is 2°. Scans (b)-(c) have the same 0.8 ^m size. In (b) the colour-scale varies over 
4nm. 

flake), we introduced a tear and fold into the sample F2 as seen in Fig. [sj^b). We see 
first that the surface under the flake is indeed free from PMMA droplets and therefore the 
topography of the flake is only influenced by the silica roughness. (This conclusion was 
also confirmed by similar measurements on other flakes). Having a flake fold allows us 
to détermine better the thickness of the flake, by measuring the step height between two 
graphene areas (as opposed to measurements of the step height between silica and graphene 
which always give a larger value of the step, ~ Inm). We find that the thickness of the 
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flake is < 0.5 nm, which confirais that the fiake is a monolayer (supporting the results of 
the quantum Hall measurements discussed in the main text). An interesting resuit from the 
tear is that the graphene fiake has a tendency to form larger ripples when detached from the 
silica surface, with a ripple height ~ 0.5 nm and width 20 nm. (The roughness of the fiake 
on the substrate is ~ 0.3 nm, see the main text.) When comparing the surface roughness of 
silica and graphene, Fig. |3](c,d), we see that the surface height variation on the clean silica 
surface is ~ 60% larger than on the graphene, i.e. graphene significantly smoothes out the 
substrate roughness. 
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